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ABSTRACT 

Based on our recent work on tidal tails of star clusters (jKiipper at al.ll201Cl[ ) we in- 
vestigate star clusters of a few 10'* M0 by means of velocity dispersion profiles and 
"^ ■ surface density profiles. We use a comprehensive set of A^-body computations of star 

clusters on various orbits within a realistic tidal field to study the evolution of these 
O ■ profiles with time, and ongoing cluster dissolution. 

f-^ ' From the velocity dispersion profiles we find that the population of potential escapers, 

C/3 I i.e. energetically unbound stars inside the Jacobi radius, dominates clusters at radii 

^ ^- above about 50% of the Jacobi radius. Beyond 70% of the Jacobi radius nearly all 

stars are energetically unbound. The velocity dispersion therefore significantly devi- 
ates from the predictions of simple equilibrium models in this regime. We furthermore 
^ ' argue that for this reason this part of a cluster cannot be used to detect a dark matter 

f~~« I halo or deviations from Newtonian gravity. 

t^^ . By fitting templates to the about lO'* computed surface density profiles we estimate 

^^D ' the accuracy which can be achieved in reconstructing the Jacobi radius of a cluster in 

'^ I this way. We find that the template of King (1962) works well for extended clusters on 

If^ . nearly circular orbits, but shows significant flaws in the case of eccentric cluster orbits. 

f^ ' This we fix by extending this template with 3 more free parameters. Our template 

^^ , can reconstruct the tidal radius over all fitted ranges with an accuracy of about 10%, 

and is especially useful in the case of cluster data with a wide radial coverage and 
for clusters showing significant extra-tidal stellar populations. No other template that 
we have tried can yield comparable results over this range of cluster conditions. All 
^^ ' templates fail to reconstruct tidal parameters of concentrated clusters, however. 

H I Moreover, we find that the bulk of a cluster adjusts to the mean tidal field which it 

experiences and not to the tidal field at perigalacticon as has often been assumed in 
other investigations, i.e. a fitted tidal radius is a cluster's time average mean tidal 
radius and not its perigalactic one. 

Furthermore, we study the tidal debris in the vicinity of the clusters and find it to 
be well represented by a power-law with a slope of -4 to -5. This steep slope we as- 
cribe to the epicyclic motion of escaped stars in the tidal tails. Star clusters close to 
apogalacticon show a significantly shallower slope of up to -1, however. We suggest 
that clusters at apogalacticon can be identified by measuring this slope. 

Key words: galaxies: kinematics and dynamics - galaxies: star clusters - methods: 
analytical - methods: A^-body simulations 



1 INTRODUCTION 

Velocity dispersion profiles and surface density profiles are 

among the most basic tools for investigating the structure 

* E-mail: akuepper@astro.uni-bonn.de (AHWK); of star clusters. However, such investigations indicate 

pavel@astro.uni-bonn.de (PK); h.baumgardt@uq.edu.au (HB); that the region around the tidal radius, at which the 

d.c.heggie@ed.ac.uk (DCH) internal acceleration of a star cluster is similar to the tidal 



2 A.H.W. Kiipper, P. Kroupa, H. Baumgardt and B.C. Heggie 



acceleration due to the galactic tidal field, is particularly 
poorly understood. 

Velocity dispersion profiles sometimes show pecu- 
liarities whic h hav e been discussed in the literature. 
iDrukier et al.l ([1993), for example, observed a flattening 
in the outer parts of the velocity dispersion profile of the 
Galactic globular cluster M15 which they interpreted as 
an effect of tidal heating by the general Gal a ctic ti de or 
by tidal shocks. IScarpa. Marconi fc Gilmozzj l|2003l ) also 
found a flattening of the velocity dispersion profile for lj 
Cen and more recently for other Galac tic globular cluster s 
like NGC6171, NGC7099 and NGC288 (|Scarpa et al.ll2007l ). 
The deviation from an expected Keplerian fall-off in the 
velocity dispersion profile occurred in all clusters at radii 
where the internal gravitational acceleration is about 
ao = 1.2 X 10~*cm s~^ and was therefore interpreted by 
Scarpa et al. a s a hint of Modified Newtonian Dynamics 
IjMilgroml Il983l ) in globular clusters. Alternatively, they 
briefly discussed the possible effect of tidal heating or a 
dark matter halo on the cluster stars. 



On the contrary, McLaughli n fc MevlanI (|2003bh found 
that by fitting a IWilsoni (jl975j) model, which has a less 
sharp cut-off at the tidal radius than the commonly used 
iKind () 19661 ) model, to the re-analyzed uj Cen data, its 
velocity dispersion profile could be explained without mod- 
ifying Newtonian gravity and without adding d ark matter. 
Similar inve stiga tions of cluster profiles , e.g. iLane et al.l 
(|2009l . l20ld ) and iBaumgardt et all (|201(]| ). were also not in 
favour of MOND. 



But it is not only the velocity dispersion profiles of 
star clusters which behave strangely at the tidal boundary; 
their surface density profiles also show peculiarities and 
sometimes controversial behaviour. 



iKind ([1963) showed that the surface density profiles of 
many globular clusters can be fitted by a simple analytical 
formula having a sharp cut-off radius, which could be 
interpreted as the tidal radius of the cluster. Later he 
derived a set of physically motivated models, with a similar 
cut-off radius corresponding to an energy cut-off in the 
energy distribution function of the cluster stars , whi ch 
provided an even better fit to cluster profiles (King|l96g). 



In contrast to that, Elson. Fall fc F reemanI (|l987h . 
and more recently iGouliermis et al. I ((2010), found that 
young massive clusters in the LMC show an exponential 
surface density profile without any tidal truncation at 
the expected tidal radius. They interpreted these find- 
ings as being due to tidal debris which was expelled at 
birth from the clusters and has not had time to disperse yet. 



ICote et al.1 



Furthermore, ICote et al.1 (120021) . 

ICarraro. Zinn fc Moni BidinI (|2007l ) and ICarrard (|2009l ') 
find the outer halo Milky-Way globular clusters Palomar 
13, Whiting 1 and AM 4, respectively, to have a significant 
excess of stars at the tidal boundary, which makes any 
fit to the surface density data very inaccurate. For all 
clusters they find the radial surface density profile, S(_R), 
to be well represented by a power-law R^ with slopes of 



about ri ~ —1.8. This excess of stars is interpreted by the 
authors as heavy mass loss in a final stage of dissolution. 
The same was found for Palomar 5, which is a well studied 
MW globular cluster clos e to the apogalacticon of its orbit 
(|Odenkirchen et al.|[20o3 ). 



Moreover, iMcLaughlin fc van der Marell ((20051) showed 
that most globular clusters of the Milky Way, the LMC and 
SMC, as well as of the Fornax dwarf spheroidal, are more 
extended than can be explained by iKing (119661^ models , 
and therefore are better represented by IWilsonI (|l975r ) 
models. In this context they emphasised the lack of physical 
understanding of this phenomenon. 



On the contrary, iBarmbv et al.l (|2009l ') found that of 23 
young massive clusters in M31 most were better fitted by 
King ( 196(j) models as these clusters do not show extended 
haloes. 

As a co n sequen ce of this lack of understanding, 
iBarmbv et al.l (i2009) asked in their investigation how 
robust the physical parameters are which were derived in 
such analyses. This they tried to estimate by analyzing 
artificial clusters in the same way as the real observa- 
tions. A simila r ana lysis has also been performed by 
iBonatto fc Bica ' (2008"), although they tested the analytic 
formula of [King, (Ji)62) with artificial observations under 
various limiting conditions. Both investigations came 
to the conclusion that physical parameters can in princi- 
ple be well recovered from such idealized mock observations. 

But probably the idealized nature of these investiga- 
tions is misleading, since both tests were performed with 
dynamically unevolved clusters. However, a self-consistent 
test of deriving physical parameters from a set of numerical 
computations of star clusters with a range of initial param- 
eters has not been performed yet. Just a few investigations 
have touched this topic so far by means of numerical com- 
putati ons (e.g. Capuzzo Dolcetta. Di Matteo, fc Miocchi 
i2005| , iDrukier et al.l |2007| . iTrenti. Vesperini fc Pasquatc 
I2OI0I ). This is due to the fact that fast codes for globular 
cluster integration like Fokker-Planck or Monte-Carlo codes 
are not able to address this problem properly, as they 
cannot follow the evolution of the tidal debris and are 
restricted to cut-off criteria in energy or angular momentum 
space. Moreover, limits in computational power prevented 
us from carrying out such investigations by means of 
coUisional A''-body codes. But the recent improvements in 
computational speed of A'^-body codes and the availability 
of accelerator hardware (GRAPEs, CPUs) now allows us to 
study the dynamical evolution of star clusters with masses 
of up to several times 10* M© . 

Thus, by computing various star clusters in a range of 
tidal conditions over several Gyr we investigate the structure 
of star clusters, and in particular the evolution of the region 
around the tidal radius - the transition region between the 
star cluster and the tidal tails - in terms of velocity disper- 
sion and surface density profiles. The paper is organized as 
follows: first, a brief introduction to the topic of potential 
escapers is given in Sec. [2] as these stars play a key role in 
this investigation. Then a few methodological remarks will 



Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters 3 



be made in Sec. [3] before we come to the velocity dispersion 
profiles in Sec. 14. fl and then to the surface density profiles 
in Sec. 14.21 In the last section we will give a summary with 
a brief discussion. 



2 POTENTIAL ESCAPERS 

A major problem in producing velocity dispersion and 
surface density profiles from observations is how to disen- 
tangle cluster members from background/foreground stars. 
By only looking at stars which are inside a projected, 
estimated Jacobi radius with similar radial velocities and 
colours it is at least possible to distinguish between the 
cluster population and field stars. Even here, however, 
the question arises whether the choice of the cut in radial 
velocity significantly influences the re sults, especially of 
the velocity dispersion profile (see iKiipper fc Kroupal 
l20ld ). Moreover, one has to be careful with how the 
Jacobi radius was estimated, i.e. was it evaluated using a 
mass estimate and was the mass estimated through the 
velocity dispersion? Or was the Jacobi radius estimated 
by a cut-off radius in the profile and is the assumption 
that this cut-off radius is equal to the Jacobi radius 
reasonable? What if the cluster is on an eccentric orbit 
about the galaxy; how does this influence the cut-off radius? 

Besides, the observer does not know if the stars in 
the sample which was extracted in this way are actually 
members of the cluster or if they are already evaporated 
from the cluster and are now part of its tidal debris, and 
just lie in projection within the cluster. 

But even if the observer could distinguish between stars 
inside the Jacobi radius and those that are beyond this 
radius, it would still be unclear if the stars in the sample are 
bound to the cluster or have already gained enough energy 
to leave the cluster but just haven't done so yet. These, 
so-called, potential escapers will have a significant infiuence 
on both the surface density and the velocity dispersion 
profile, if there is a non-negligible fraction of them in the 
cluster, as these stars will make the profiles deviate from 
any theory which does not take them into account - which 
most theoretical approaches do not. 

In numerical modelling of star clusters we have detailed 
phase space information on every single star in the compu- 
tation. Mainly by numerical investigations it has been found 
that the escape process through which stars escape from a 
cluster in a constant tidal field, e.g. on a circular orbit about 
a galaxy, is divided into two steps. First the stars get un- 
bound via two-body relaxation, thus on a relaxation time 
scale (ignoring constants). 



f-rei O^ -'• t^cr^ 



(1) 



where A^ is the number of stars in the cluster and tc- is 
the crossing time of the cluster (note that we neglected the 
slowly varyi ng Coulomb logarith m in eq. [U for further de- 
tails see e.g. iHeggie fc Hu1il2003l 'l. In the second step, these 
energetically unbound stars, or potential escapers, with spe- 
cific energy E higher than some critical escape energy Ecru, 



escape fr om the cluster on an esca pe time scale which is 
given by (JFukushige fc Heggiell200d . equation 9 therein) 



T-esc CX Lev 



( E — Ecrit \ 
V Ecrit ) 



(2) 



where Ecrit is the critical energy at the Lagrange points, 
i.e. at the Jacobi radius. From this it follows that the excess 
energy can be written as 

E — Ecrit 



0^ \yc.r j^esc) 



1/2 



On the other hand, the relation 



Ot ytesc /i^rel ) 



1/2 



E — Ecr 
^crit 

holds, whence we find that 

'^esc ^ \''cr '^rel ) 



(3) 



(4) 



(5) 



Hence, a fraction of stars with excess energy {E — 
Ecrit) /Ecrit escapes on a time scale tcsc and the time scale 
of mass loss, e.g. the time on which a cluster dissolves, tdiss, 
is given by 

^E-Ecr^Y^ (6) 



td: 



OC tc 



Ecrit 






(7) 

OC trcl N-^'^ (8) 

which was also found bv iBaumgardtl (|200l|). This result is 
in contrast to the 'classic' picture where the dissolution 
times of clusters scale w ith the relaxation time (see e.g. 
iBinnev fc TremaindlJOO^ ), i.e. 



f^diss CX Lrel 5 



(9) 



and demonstrates the importance of potential escapers 
on the dissolution process of star clusters. The remaining 
question is whether the population of potential escapers is 
large enough that it also has an infiuence on the velocity 
dispersion and surface density pro files. For clus t ers of a few 
10'' stars in a constant tidal field IBaumgardtl l|200ll ) finds 
about 10-20% of all stars wi t hin th e Jacobi radius to be 
potential escapers, Ijust et al.l (|2009l ) even find the fraction 
of potential escapers for sim ilar clusters to be 1/3 of the 
cluster population. Moreover. IBaumgardtl l|200ll ) found that 
for clusters in a constant tidal field the fraction of potential 
escapers varies with the number of stars within a cluster 
approximately as N~^''^. Thus, this fraction should still be 
significant for high- A'' globular clusters. 

Star clusters in time-dependent tidal fields have not 
been investigated in this respect yet, but as tidal pertur- 
bati ons tend to increase the energ y of the stars in a cluster 
(e.g. Gne din. Lee fc Ostrikerll9 99!) , the population of poten- 
tial escapers should be even larger in such clusters. Thus, 
potential escapers should have a significant effect on both 
the velocity dispersion and the surface density profiles of all 
kinds of star clusters. 



3 METHOD 

When analyzing velocity dispersion and surface density 
data of star clusters, confusion over nomenclature often 
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arises. Therefore we strictly stick to the foUowing three 
terms: profiles, models and templates. The first will be 
used for the data, which in our case comes from A^-body 
computations but might also originate from observations. 
The term models will be used for physically moti vated 
distrib utio n funct i ons s uch as those proposed by iKind 
l| 19661 ') or IWilsonI (llOTSl ). The word templates covers the 
analytic expressions which have been found empirically and 
which do not originate from a physical derivation by being 
solutions of the Collisionless Boltzmann Equation (e.i 
iKindl 19621 : lElson. Fall fc Freeman|[l987l : iLauer et al.lll995 
The word template hereby emphasises their artificial nature. 

Further confusion arises about the term tidal radius. 
Again there is a physically motivated and an empirical 
version. The former is the radius at which the internal 
gravitational acceleration equals the tidal acceleratio n from 
the host galaxy (see e.g. iBinnev fc Tremaind I2OO81 I. This 
will be named Jacobi radius in our text. The empirical tidal 
radius is defined through the cut-off radius of the models 
and templates which are fitted to the profiles. We will refer 
to this cut-off radius as edge radius here and would like 
to emphasise that Jacobi radius and edge radius are two 
different concepts and therefore are not necessarily equal. 
In fact, we will show that the two are often significantly 
different and that conclusions from a fitted edge radius 
should be drawn with caution. 

Note also that the results of this investigation cannot 
be easily scaled to more massive globular clusters, as the 
relative importance of tidal features depends on the mass of 
the cluster, since the mass in the tidal debris is proportional 
to the mass-loss rate of the cluster and the mass-loss rate 
does not sc ale linearly with cluster m ass, M, but rather 
with M^/" JBaumgardt fc Makindl2003l '). 

Moreover, even though we have some concentrated 
clusters in our sample this investigation focuses on less 
concentrated clusters, with a ratio of projected half mass 
radius to Jacobi radius betwee n 0.1 < Rhp/Rj < 0.3 
because iBaumgardt et al.l ([20101) showed that clusters with 
ratios of Rh v/Rj < 0.05 cannot be pr operly fitted by either 
iKing l| 19621 ) templates or iKina (|l966l ) models as their ratio 
of Rhp/Rj is too high. Therefore, our conclusion s mainly 
hold for less concentrated clusters which Baumga rdt et al.l 
l|201Cl) named the extended cluster population, and which 
are likely to be on the main sequenc e of cluster evolution 
(|Kiipper. Kroupa fc Baumgardt! 120081 ) , the final and uni- 
versal stage of cluster evolution after core collapse, where 
the ratio Rhp/R.j only depends on the mass of the cluster. 

In the following sections we will first give some details 
about our data set before we explain how we extract the 
velocity dispersion and surface density profiles. Thereafter 
we briefiy describe the analytical templates which are used 
in this investigation and specify our fitting method. 



3.1 Data set 

The velocity dispersion and surface density profiles are 
mainly taken from o ur recent investigat ion of star clusters 
and their tidal tails IjKiipper et al.ll2010l ). For this analysis 



Table 1. Overview of all computed models. R?f° gives the 
apogalactic radius of the cluster orbit and e the corresponding 
eccentricity. The latter is defined as e = {R'q° — FSq^^) / {R'q° + 
RFq") with R'q° being the apogalactic radius of the cluster orbit, 
and R'^^' being the perigalactic radius. Rf^p/ Rj gives the initial 
concentration of the model, whereas R}^ is the initial half-mass 
radius. The clusters are divided into three categories: extended 
clusters in a constant tidal field (AO-DO), extended clusters in 
time-dependent tidal fields (A1-D3) and compact clusters (AOc- 
DOc). Note that each model has initially N = 64fc stars, a lower 
stellar mass limit of 0.1 Mq and an upper limit of 1.2 Mq. 



Name 


R''^" [kpc] 


e 




Rhp 
Rj 


Rh [pc] 


AO 
BO 
CO 
DO 


4.25 
8.50 
12.75 
17.0 


0.00 
0.00 
0.00 
0.00 




0.15 
0.15 
0.15 
0.15 


5.8 
8.8 
11.7 
14.3 


A1/A2/A3 
B1/B2/B3 
C1/C2/C3 
D1/D2/D3 


4.25 
8.50 
12.75 
17.0 


0.25/0.50/0 
0.25/0.50/0 
0.25/0.50/0 
0.25/0.50/0 


75 
75 
75 
75 


0.15 
0.15 
0.15 
0.15 


5.8 
8.8 
11.7 
14.3 


AOc 
BOc 
COc 
DOc 


4.25 
8.50 
12.75 
17.0 


0.00 
0.00 
0.00 
0.00 




0.08 
0.06 
0.05 
0.04 


3.0 
3.0 
3.0 
3.0 



we computed a comprehensive set of star clusters on various 
orbits about a Milky-Way potential. 



In addition to the 16 clusters taken from lKupper et al.l 

([2010!), we computed 4 concentrated clusters which are of 
specific interest for this investigation. All clusters were set- 
up using the publicly available code McLustefij (Kiipper et 
al., in prep. ), and integrat ed over time with the TV-body code 
NBODY4 iAarsethll2003l) on the GRAPE-6A c omputers at 
AlfA Bonn (|Fukushige. Makino fc Kawaill2005l ). The code 
was modified so that the Milky- Way potential suggested by 
[Allen fc Santillan (1991) could be used. Due to its analytic 
form this potential is useful for calculating the Jacobi radius, 
Rj, of the star clusters. Through out the compu t ation s Rj 
was evaluated using equation 12 of [Kiipper et al.l ([2OI0I ). i.e. 



R.J 



GM 



0,2 - d^^/dRl ' 



(10) 



where G is the gravitational constant. A/ is the cluster 
mass, il is the angular velocity of the cluster on its orbit 
about the galaxy, and d^^/dR^ is the second derivative of 
the gravitational potential of the galaxy with respect to the 
galactocentric radius, Rq. We always determine the Jacobi 
radius iteratively, assuming first a mass of the cluster, then 
applying eq. 1101 With this estimate we determine the mass 
of the cluster again, counting all masses around the cluster 
centre within R,j, and compute the Jacobi radius once 
more, and so on until the value of R,j converges. 

The clusters were set up using a Plummer density 
distribution with a ratio of half-mass radius to Jacobi radius 
of Rh/R.j = 0.2, i.e. a ratio of projected half-mass radius 



www . astro . uni-bonn . de/'akuepper/mcluster/mcluster . html 
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to Jacobi radius of Rhp/Rj = 0.15. This choice resembles 
the extended globular clusters of the Milky Way, which are 
less concentrated, and stand in contrast to the compact 
clusters, which are deeply embedded within their Jacobi 
radius and for which ti dal influences are rather negligible 
(JBaumgardt et al.l[201ol ). 

All clusters had initially 65 536 stars wh ich were 
drawn from the canonical IMElj (JKroupal 1200 ih ranging 
from 0.1 M0 to 1.2 M©, resulting in a total initial mass 
of about 2OOOOM0. All clusters were modelled without 
stellar evolution and without primordial binaries. They 
were computed for a time span of 4 Gyr, if they did not 
dissolve before reaching this age. An overview of all clusters 
in our sample is given in Tab. [T] 

For all clusters we produced two-dimensional snapshots 
of the stars projected onto the orbital plane every 10 Myr; 
thus, as each of the 20 clusters is modelled over 4 Gyr, we 
have about 20 x 400 ~ lO** representations of star clusters. 
In each snapshot we deteriuine the centre of the cluster stars 
within the Jacobi radius following Cascrtano & Hut (1985) 
and bin the stars around this centre in 50 annuli of equal 
logarithmic width between 1.0 pc and 500 pc to produce 
the velocity dispersion and surface density profiles. By go- 
ing out to 500 pc we are able to show not only the region 
close to the Jacobi radius but also a considerable part of the 
tidal tails, which have not been investigated in this respect 
yet, but show an interesting, variable behaviour, depending 
on the cluster orbit, which can even influence the velocity 
dispersion and surf ace density within the Jacobi radius (see 
iKiipper et al.ll2010r ). For most of the investigation, however, 
we concentrate on the inner 100 pc. 



we will investigate the effect of foreground/background 
stars, i.e. stars in the tidal debris, on the observed velocity 
dispersion. Finally we measure a taking only the bound 
stars into account to test if energetically unbound stars have 
a significant influence on the observed velocity dispersion. 

For this purpose we have to search for potential escapers 
and remove them from the sample. Thus, bound stars are 
deflned as stars which are inside the Jacobi radius and which 
have a Jacobi energy smaller than 



^crit — 



3GM 
' 2Rj ' 



(13) 



which is the critical energy at the Lagrange points, i.e. at the 
Jacobi radius. We therefore compute the Jacobi energy, E, 
of each star with mass m in a co-rotating reference frame, ac- 
co rding to the n ear-field approximation which can be found 
in ISpitzerl(| 19871 ): 



^ mv GMm m^o / 



' 3a:: 



(14) 



where the first term is the kinetic energy of the star, 
the second term is its energy in the gravitational field 
of the cluster with mass M at the star's radius r, and 
the third term is a combination of the centrifugal and 
tidal potentials. In this last term z gives the distance 
of the star from the cluster centre perpendicular to the 
cluster's orbital plane in the galaxy, whereas x denotes 
the coordinate of the star relative to the cluster centre, 
along the axis pointin g to the galactic centre (for a detailed 
description see e.g. iFukushige fc Heggid l2000h . We are 
well aware of the fact that this approximation holds only 
for nearly circular cluster orbits in a point-mass galac- 
tic field, and take it only as a first-order approximation here. 



3.2 Velocity dispersion profiles 

In each annulus the velocity dispersion along the line-of-sight 
(los), a, is determined using 



iV 
O-^ = — ^ («i - vf = W2 



(11) 



where Vi is the velocity of the i-th star along the line-of-sight, 
V is the mean los velocity of the A'^ stars in the annulus, 
and w^ is their mean squared los velocity. In addition, we 
estimate an uncertainty for each velocity dispersion. 



1=1 



[{vi - vf - cr^) , 



(12) 



which originates from a Taylor expansion of the standard 
deviation of the velocity dispersion. 

The velocity dispersion is first measured by taking 
into account all stars in the snapshot, independent of their 
projected distance from the cluster centre, which is what 
an ideal observer would see. Then it is measured for only 
the stars within the Jacobi radius, i.e. r ^ R,j. In this way 

^ The canonical IMF has a slope of oi = 1.3 for stellar masses 
m = 0.08 — 0.5 Mq, and the Salpeter slope a2 = 2.3 for m > 
O.5M0. 



Note that we remove binaries from our data and 
replace them by their centre-of-mass particles, since 
we want to focus on the effect of potential escapers 
on the velocity dispersion. In observational data, bina- 
ries, of course, c an play a major role, as has recently 



been shown by Gieles. San a & Portegies Zwart (200! 



iKouwenhoven fc de GriisI (12009) and Kiippcr & Krou|3; 



( 20101 ) but only for clusters with a velocity dispersion of 
less than about 5-10 km/s. 

3.3 Surface density profiles 

In each annulus the surface density, E(_R), is determined by 
counting the number of stars. A'', in the given (projected) 
radial range [R — AR, R] , and dividing by the surface area, 
A, of the given annulus, i.e. 

^(^) - f = 2niR^-^R-ARry (15) 

Moreover, for each annulus an uncertainty is estimated by 
the square root of the number of stars in the annulus divided 
by its surface area. 



AE(i?) = 



2t: {R2 ~ {R ~ AR)'^) 



(16) 



First the surface density profiles are evaluated by taking into 
account all stars in each sample, and then additionally only 
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Table 2. Overview of the templates used in this investigation. 
nf gives the number of free parameters, where the +1 denotes an 
additionally assumed constant background. As each template was 
originally designed with the focus on specific components of star 
clusters, this table shows for each template which of the three 
basic components it can handle: core, bulk and tidal debris. 



Name 



bulk 



tidal debris 



King 


3+1 


flat 


separated 


n/a 


EFF 


3+1 


flat 


not Sep. 


power-law 


Nuker 


5+1 


cuspy 


not Sep. 


power-law 


gen. Nuker 


8 


cuspy 


separated 


power-law 


KKBH 


6+1 


cuspy 


separated 


power-law 



10' 



KKBH (7=0.0) 

KKBH (7 =0.2) 

f1(R) 

f2(R) 

b 




1000 



R[pc] 



Figure 1. Sketch of the KKBH template for two different val- 
ues of the core slope 7. Also shown are the two sub-components 
for the bulk (fl{R), eq. I22I I. which resembles the King (1966) 
te mplate, and the tidal debris (/ 2(fl), eq. I23I I. which resembles 
the lElson. Fall fc FVeemar] l ll987t) template. The break radius be- 
tween the two components is at fiRt, where fi is found to be 0.5 
in most cases, i.e. for radii larger than 50% of the edge radius the 
surface density profiles of the clusters follow a power-law. More 
flexibility for large radii is added to most templates through a 
constant background b. 



the bound stars, to show where and in which way potential 
escapers influence the profile. 



3.4 Analytical templates 

Ev en thoug h physically motivated models like iKind (|l966l ) 
or IWilsorJ (|l975. ) have often proved to more generally 
represent observ ations of globular clusters than an alytical 
templates (e.g. iMcLaughlin fc van der Marell |2005| ). they 
have one decisive flaw: there is no simple analytical de- 
scription which can be quickly fitted to star cluster profiles. 
Thus, due to their simplicity, analytical templates are still 
in use, especially for larger data sets of extra-Galactic 
clusters. 

Moreover, analytical templates have the advantage 



that they can be easily set up and extended, according 
to the needs of the corresponding investigation. We found 
that star clusters consist basically of three parts: core, bulk 
and tidal debris. The templates which have been set up 
in the past were adjusted according to the differen t foci 
of th e corresponding investigations. For example, iKind 
l|l962r ) was mainly interested in the bulk of the cluster, 
and so the templ ate mainly consists of a bulk and just a 
flat core, whereas lElson. Fall fc FreemanI l|l987i ) focused on 
the tidal debris of young clusters, and so their template 
has a flat core, no separated, explicitly defined bulk but a 
power-law distribution of tidal debris. Therefore, we will 
study some existing templates in a systematic way to show 
their advantages and their flaws, and finally design a more 
general template. An overview of all templates is given in 
Tab. [3 

The 'classic' templates we use, their abbreviations for 
the rest of the text, and their fit parameters are the follow- 
ing: 

• King denot e s the empirical fitting formula which was 
devised by iKingj (|l962l ) and which works surprisingly well 
for many extended globular clusters (i.e. Rhp/Rj > 0.05) 
and open clusters. It has the form 



f{R)=k 



v^l + {R/Rcf x/l + {Rt/Rc 



+ b (17) 



for R < Rt and f{R) = b for R ^ Rt, where fc is a scale 
factor. Re is a core radius, and Rt is often denoted as the 
tidal radius - in our nomenclature it is the edge radius, the 
limiting radius between the bulk of a cluster and its tidal 
debris. The King template therefore consists of a flat core 
and a bulk, but has no term for the tidal debris. Note that 
we allow for a constant background b instead, as would 
be done in real observations even though we have no real 
background but only the tidal debris of the clusters. We 
found this constant to significantly improve the quality and 
stability of the fit. 

• EFF is the empi ri cal te mplate used in the work of 
lElson. Fall fc Freemaal (|l987l ) on young LMC clusters of 
about the same mass as the clusters in our sample. They 
found that the outer parts of those clusters were poorly fit- 
ted by the King template and adopted this template which 
has no edge radius but falls off like a power-law: 



f{R) = fc 



-(I 



-,,/2 



(18) 



where fc is a scale factor. Re is a core radius, & is a constant 
background and 77 is the slope of the template for radii much 
larger than the core radius. The EFF template was designed 
to fit observations of clusters with pronounced tidal debris, 
and therefore consists of a flat core and an unseparated bulk 
and tidal debris which are represented by a single power-law. 

• Nuker denotes the template which was adopted by 
iLauer et al.l (|l995l ) for elliptical galaxies and which has also 
occasionally been used for star clusters. It is also a power- 
law but shows more flexibility than EFF and can also flt a 
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cluster with a non-flat core, e.g. a core-collapsed cluster or 
a cluster with a massive central black hole. It has the form 



fiR) = k2- 



1 + 



\Rc) 



+ b, 



(19) 



where fc is a scale factor, Re is a core radius and 6 is a 
constant background. 7 gives the power-law slope inside 
Re, whereas rj gives the slope for radii larger than the 
break radius Re- With the factor a, the smoothness of the 
transition between the two slopes is defined. Thus, the 
Nuker template consists of a fiexible core and a power-law 
tidal debris. 

Additionally, we set up two templates which have a fiex- 
ible core plus an edge radius, i.e. a well-defined bulk, but can 
also fit the region beyond the edge radius in a more sophis- 
ticated way than by just a constant background. 

• Generalised Nuker, gen. Nuker, is the same as Nuker 
but allows for a different power-law slope for radii larger 
than a s econd break radius. The tem plate has been intro- 
duced bv lvan der Marel &: AndersonI (2010) to get a smooth 
representation of the surface density profile of uj Cen. It has 
the form 



fiR) 



'-'Itj 


-7 - 
1 


h(l)l 


s 



R\°' 



(f) 



(20) 



where the parameters are the same as in Nuker, though 
without the constant background, but instead with an 
additional second break radius, Rt, a slope for radii larger 
than this radius, e, and a factor 5 which determines the 
smoothness of the tran sition between e and rj. I n their 
investigation of uj Cen, Ivan der Marel fc AndersonI (|2010|) 
found that the two break radii in this template roughly 
correspond to the core and Jacobi radius of the cluster. In 
our nomenclature their template consists of a fiexible core, 
a bulk and a tidal debris, which are separated by two break 
radii: the core radius and the edge radius. 

• Finally, we set up our own template which is tailored 
to the purpose of this investigation. KKBH is based on the 
King template but modified in two steps: first we modified 
the King template in much the same way as the Nuker tem- 
plate is modified with respect to EFF, i.e. we allow it to 
have a power-law cusp in the core, but without changing 
the behaviour in the outer parts of the template: 



f{R) = k 




+ b (21) 



for R < Rt and f{R) — b for R ^ Rt, with the same 
parameters as King plus an additional core power-law slope 
7. For 7 = it gives the original King template. This 
template therefore consists of a fiexible core and a bulk but 
still has no tidal debris term, yet. 

For this reason we modified it further with an additional 



extra-tidal component in the form of a power-law slope. For 
the core and the inner part of the bulk of the cluster it 
behaves like the previous template, but at a radius which is 
a fraction ^ of the edge radius the template changes abruptly 
into a power-law, and thus behaves like EFF (see Fig. [ij . It 
is given by 




MR) 

y'l + iR/Rcf \/l + (Rt/Rc 
for radii smaller than jj,Rt, and 



f2iR)^Mf^Rt) 



1 + 



R 

fj-Rt 



-J7/64 



(22) 



(23) 



for R^ fiRt, where the parameters are the same as for the 
previous template, except for the new ones in the function 
/2(-R). The exponent 64 in f2{R) causes the template to 
change abruptly into the power-law slope of 77 at a fraction /i 
of the edge radius. With an additional constant background, 
6, the complete function looks as follows when using ternary 
operators 



f{R) = (R< fiRt 



fi{R) + b 



MR) + b) 



(24) 



(meaning: if R is smaller than fiRt use fi{R) + b, else use 
f2{R) + b). Defined in this way, the KKBH template has 
a fiexible core, a well-defined bulk and a power-law tidal 
debris (see Fig. [T]). A brief explanation on how to use the 
KKBH template with gnuplot is given in Appendix A. 

The three additional parameters of KKBH may not be 
useful in all applications but are motivated through the fol- 
lowing facts: 

- Mass segregation or the presence of an intermediate- 
mass black hole can have a significant infiuence on the 
slope of the core profile. This cannot be fitted and thus 
quantified by the original King template, and therefore we 
add the parameter 7 to the template. 

- In most applications the star counts in the outer 
parts of surface density profiles drop quickly, while 
their uncertainties grow, so that a template, and hence 
the edge radius, is mo stly fitted by the inner profile 
JBaumgardt et al.l I2OICI ). Hence, by allowing a smooth 
template for the inner profile and an independent power- 
law for the outer profile we decouple the two parts and 
look for the transition point between the two, which is set 
by the second additional parameter fi. 

- Potential escapers and background/foreground stars 
as well as the tidal debris of the cluster i nfiuence the tem- 
plate fits in the outer parts of profiles (jC otc ot al .1 120021 : 
ICarraro. Zinn fc Moni BidinI l2007l : ICarraro 2009 ') . Thus, 
the tidal debris is not reflected properly by templates like 
King, and so the results of template fltting are not inde- 
pendent of the presence of debris. Furthermore, by mea- 
suring the slope of the surface brightness profile outside 
the edge radius we aim to be able to deduce information 
on the orbital phase of the cluster, and therefore we add 
the parameter 77. 
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3.5 Fitting method 

We fit the above templates to the data by finding a minimum 
in X using a nonlinear least-squares Marquardt-Levenberg 
algorithm. Here we define x^ as a measure of the differences 
between the template, f{R), and the data, E(J?), weighted 
by the uncertainties in the surface density data, AE(_R), 
summed up over all data points, i.e. 



X 



E 






(25) 



To compare the fits to each other and to account for the 
different complexity of the templates we calculate a reduced 
X^ which is given by the above x^ divided by the number of 
degrees of freedom, n. 



(26) 



Area 5 

n 

where n is given by the number of data points, Ud, minus 
the number of parameters of the given template, n/ , minus 
1 (as one parameter can be fixed by the mean value of the 
data points), i.e. 

n = Ud — rif — 1. (27) 

The resulting values of Xred should ideally be close to one, 
as too low a value would imply that the template is too 
complex for the underlying data. 



Each of the clusters in a constant tidal field has a 
steadily declining velocity dispersion within the Jacobi 
radius (green/red), although the clusters in the stronger 
tidal fields already show an onset of a flattening of the 
profile at the Jacobi radius. Beyond the Jacobi radius 
the velocity dispersion fluctuates strongly as the number 
statistics in these outer bins is quite low and a single 
fast escaper may dominate the velocity dispersion in the 
corresponding bin (red). 

From the differently coloured curves we can see that 
the stars which lie within the projected Jacobi radius, 
but actually are in the foreground or background of the 
cluster within its tidal debris, hardly influence the velocity 
dispersion of the clusters (green). Only in some of the 
graphs can one make out a small contribution in the very 
outermost bins; everywhere else the green data points lie 
on top of the red ones. 

By contrast, potential escapers have a significant infiu- 
ence on the profiles (green data vs. blue data). By definition 
the velocity dispersion of the bound stars approaches zero 
at the Jacobi radius, where the tidal forces are equal to the 
gravitational acceleration of the cluster. But we see a signif- 
icant infiuence already at about 50% of the Jacobi radius. 



Each template is fitted to each cluster for five differ- 
ent ranges of radii to see if the goodness of the fit changes 
for different coverages of the cluster and its debris, which 
also implies a differing number of data points. Therefore we 
choose fixed radial intervals of 1-25 pc, 1-50 pc, 1-100 pc, 
1-200 pc and 1-500 pc. Since the number of data points be- 
tween 1-500 pc is fixed to 50 for all clusters, we get [25, 31, 
37, 42, 50] data points in the given radial ranges. As the 
clusters have the same initial mass but orbit the galaxy at 
different galactocentric radii, these five intervals will cover 
quite different parts of the clusters. For instance a cluster 
of 2OOOOM0 has a Jacobi radius of about [27, 40, 52, 63] 
pc for a circular orbit at [4.25, 8.5, 12.75, 17] kpc (eg. fTO)) . 
respectively. 



4 RESULTS 

4.1 Velocity dispersion profiles 

4-1.1 Constant tidal fields 

A typical sample of velocity dispersion profiles is shown in 
Fig. [2I In each graph the velocity dispersion is plotted i) for 
all stars in the sample (red), ii) for all stars that lie within 
the Jacobi radius (green) , and iii) for all bound stars (blue) . 

The first column shows the clusters in a constant tidal 
field, i.e. on circular orbits. Comparing the graphs for the 
different galactocentric distances shows that, as expected, 
the clusters at larger galactocentric distances are more 
extended. Since they all have the same mass of about 
15 000 Mq in the snapshots, the central velocity dispersion 
decreases with Re, whereas the size of the ffat part within 
the core radius increases. 



4-1.2 Time- dependent tidal fields 

If we take a look at columns 2 and 3 of Fig. [2] we can see 
the effect of time-dependent tidal fields on the velocity 
dispersion profiles. With increasing eccentricity the extent 
and onset of the population of potential escapers stays the 
same (~ 0.5-R.7), since in the snapshots all clusters are at 
apogalacticon. But from the green curve we can see that 
the energy of the potential escapers increases, which leads 
to a fiattening of the veloci ty dispersion profiles. Th i s is in 
agreement with findings by ICnedin. Lee fc Ostrikeij (119991 ) 
who showed that stars at large radii are most affected by 
tidal perturbations like disk shocks or pericentre passages 
and that stars, on average, gain energy throug h such events. 
Furth ermore, it supports the assumption of iKiipper et al.l 
l|201Cl) that stars do not escape from the cluster immediately 
after a tidal perturbation, but can remain in the cluster for 
several cluster orbits about the galactic centre. 

In Fig. [3] we show a time series of a cluster on an orbit 
with an eccentricity of 0.75 and an apogalactic radius of 17 
kpc. The time series starts 210 Myr after the beginning of 
the computations, where the cluster is just about to finish 
its first revolution about the galactic centre and has a mass 
of 16 2OOM0. At 310 Myr (middle panel) the cluster is 
at perigalacticon with a mass inside the Jacobi radius of 
7 700Mq, and at 410 Myr it is again at apogalacticon with 
13 500Mq left. 

We see that for the largest part of the orbit, from 250 
Myr to 370 Myr {Re < 15 kpc) the velocity dispersion pro- 
file appears to be more or less fiat. For the perigalactic part 
of the orbit (300 Myr - 320 Myr, Rg < 5 kpc) a good fraction 
of the cluster stars even lies beyond the Jacobi radius (about 
50% of the apogalactic cluster mass). Nevertheless, even af- 
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Figure 2. Velocity dispersion profiles of 12 different clusters (AO-2, BO-2, CO-2, DO-2). The apogalactic distance and eccentricity of the 
clusters' orbits are given in the graphs. At the given snapshot each cluster is at apogalacticon and has a mass of about 15 000Mq. The 
corresponding Jacobi radius can be deduced by the outermost green data point. Two effects are visible: first, for increasing galactocentric 
distance the extension of the cluster, and hence also the flat part inside the core radius, both grow as the Jacobi radius grows; this results 
from setting up the cluster half-mass radii to scale with Rj (Sec.|3TTJ- This also implies a lower central velocity dispersion. Second and 
more important, for increasing eccentricity the effect of potential escapers (green points vs. blue points) on the flatness of the velocity 
dispersion profile outside the core radius grows. Corresponding surface density profiles are given in Figs.[5]&|9] 
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Figure 3. Velocity dispersion profiles for the cluster D3, which is on an eccentric orbit with eccentricity 0.75 and an apogalactic distance 
of 17 kpc. The snapshots show the evolution of the profile with time during one revolution about the galactic centre, from apogalacticon 
to apogalacticon. In each graph the corresponding simulation time and galactocentric radius are given. In the middle graph the cluster 
is at perigalacticon. Here the Jacobi radius is at its minimum (~10 pc) and most of the cluster is unbound. For a large fraction of the 
orbit the velocity dispersion profile is approximately flat. Corresponding surface density profiles are given in Fig, ill! 



ter such a dramatic pericentre passage the cluster's velocity 
dispersion profile returns to nearly its original shape. 



near core collapse. Even so, there is no significant difference 
to the other clusters with respect to the potential escapers. 



4-1.3 Concentrated clusters 

In Fig. Uthe velocity dispersion profile for a cluster at 4.25 
kpc is shown. It had an initi al (3-dimensional) half-m ass ra- 
dius of 3 pc, and is thus what lBaumgardt et al.l l|2010r ) would 
classify as a compact cluster. At the time of the snapshot 
the cluster has 15 000 M0 and still has a ratio of projected 
half-mass radius to Jacobi radius Rhp/Rj = 0.08. The pro- 
file does not show a well-defined core in the range of the data 
points, Eis all other profiles do, but rises steeply in the cen- 
tre. This is a clear indication of the fact that the cluster is 



4-1-4 Conclusions 

From the velocity dispersion profiles studied in this section 
we can conclude that the population of potential escapers 
has a significant infiuence on the kinematical structure of a 
star cluster. The influence of this population on the velocity 
dispersion profiles grows with increasing eccentricity of 
the cluster orbit. In such cases an observer would have 
to be very cautious when deriving a dynamical mass out 
of velocity dispersion measurements as it would tend to 
be overestimated, the possible degree of overestimation 
depending on the cluster characteristics such as its orbit 
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Figure 4. Velocity dispersion profile for the concentrated cluster 
with a galactocentric radius of 4.25 kpc (AOc). At the time of the 
snapshot the cluster has a mass of about 15 000 Mq and a ratio 
of projected half-mass radius to Jacobi radius of 0.08. Compared 
to the other clusters at the same stage of evolution, the velocity 
dispersion rises steeply in the centre due to core collapse, but 
shows a similar behaviour in the outer parts. A corresponding 
surface density profile is given in Fig. 1131 



and its velocity dispersion. 

Furthermore, we saw that the velocity dispersion profile 
deviates from the expected Keplerian profile at a radius of 
about half the Jacobi radius. Since most globular clusters in 
the Milky Way experience a gravitational acceleration which 
is about the size of GMcai/ Roai — 1-2 x 10~ cm s~ ~ ao 
one may easily be misl ed to deduced dev iations from 
Newtonian gravity (e.g. IScarpa et al.l 120071 ). To reduce 
the influence of potential escapers, velocity dispersion 
measurements should be made within 50% of the Jacobi 
radius. Nevertheless, extra caution has to be taken for 
clusters on eccentric orbits. 

Moreover, we can see the general structure of star clus- 
ters already in the velocity dispersion proflles. There is a 
core which can be fiat or cuspy. Then there is a bulk of 
bound stars whose limiting radius is somewhat smaller than 
the Jacobi radius, and finally there is the tidal debris whose 
onset lies within the Jacobi radius as a consequence of po- 
tential escapers. Most of this cluster structure is temporarily 
erased during pericentre passages for eccentric cluster orbits. 
By contrast, the concentrated clusters show a different clus- 
ter structure, as they do not seem to show a distinct core 
and bulk within the range of the data points. 



25 

20 

15 

X 10 

5 





'rrr n ' i 


iP 


Nuker ri ;; 

Gen. Nuker ri ;; i 

Gen. Nuker E : \ ;l 




1 1 


1 '; 



1000 



2000 
t [Myr] 



3000 



4000 



Figure 6. Evolution of the slopes of the tidal debris of the tem- 
plates which can fit such a component, i.e. the surface density 
in the tidal debris falls off with R~^ , for the cluster on a cir- 
cular orbit at 8.5 kpc (BO). Since the gen. Nuker template has 
two power-law slopes we show both here. Note that the Nuker 
and the gen. Nuker slopes are smoothed, as the variations be- 
tween the individual time steps are very large. The other slopes 
are not smoothed since they are better behaved. The EFF and 
KKBH slopes show a rise in the beginning of the simulation which 
corresponds to an initial increased mass-loss rate due to primor- 
dial escapers. The Nuker and gen. Nuker slopes give no reliable 
information on the tidal debris slope. 
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Figure 7. Evolution of the Jacobi radius Rj of the cluster on a 
circular orbit at 8.5 kpc (BO), evaluated using eg. 1101 (red dotted 
line), and the edge radii of the fitted templates for a fit range 
of 1-100 pc. None of the edge radii reproduces the Jacobi radius 
well over the whole simulation time. King shows considerable evo- 
lution with respect to the Jacobi radius. The edge radius of the 
gen. Nuker template is smoothed since it fluctuates a lot, or some- 
times does not converge properly (see text). In cases when it con- 
verges, it mostly gives a radius which is about 50% larger than 
Rj. KKBH has about the same values as King, though reduced 
by another 10%, but is less influenced by the initial mass loss due 
to primordial escapers, i.e. is more stable. 



4.2 Surface density profiles 

4-2.1 Constant tidal fields 

As an example to demonstrate the fitting process, in Fig. [5] 
the surface density profile of the cluster on a circular 
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Figure 5. Surface density profiles of the cluster in a constant tidal field at a galactocentric distance of 8.5 kpc (BO). The upper panels of 
the first five graphs show the data and the template fits according to Sec. l3.4l for the five ranges of l-[25, 50, 100, 200, 500] pc. The lower 
panels of the first five graphs give the relative deviation of the templates from the data. Vertical dashed lines mark the fitting range. The 
theoretical Jacobi radius can be read off from the outermost green data point (bound stars). The lower right panel shows the reduced 
X^ for all of the five ranges, averaged over all snapshots of this cluster. Error bars give the standard deviation of Xred' ^ corresponding 
velocity dispersion profile is given in Fig. [21 
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Figure 8. Ratio of the fitted edge radii to the Jacobi radius av- 
eraged over the whole simulation of the cluster on a circular orbit 
at 8.5 kpc (BO) and shown for all 5 fit ranges. Error bars give 
the standard deviations from the mean. The King template re- 
produces the Jacobi radius to about 10% accuracy for sufficiently 
large fit ranges. Gen. Nuker gives values 2-3 times larger than 
the Jacobi radius and shows large fluctuations. KKBH shows a 
constant value over all fit ranges which is about 20% smaller than 
the Jacobi radius. Also shown are the edge radius of the bound 
stars, which lies at about 70% of Rj, and the break radius of 
the KKBH template at which the profile of the cluster becomes 
dominated by potential escapers. 



orbit at 8.5 kpc is shown when the cluster has a mass of 
I5OOOM0, just as was shown for the velocity dispersion 
in the 4th panel of Fig. [21 In the upper panels of the first 
five graphs the cluster data and the fits of the analytical 
templates from Sec. 13.4! are shown for the five different 
ranges of l-[25, 50, 100, 200, 500] pc. In the lower panels 
of these graphs the relative errors can be seen. The lower 
right panel gives the reduced x^ for each template for each 
fit range averaged over all 400 snapshots of the cluster. 

We see that for the smallest fit range (1-25 pc) all tem- 
plates fit about equally well, whereas for the larger ranges 
significant differences in the quality of the fits appear. This 
is due to the fact that the smallest fit range only covers the 
core of the cluster, which in this case is flat, and a part of 
the bulk. Such a profile can be fit by all templates in the 
sample. But already the second smallest fit range of 1-50 
pc covers part of the tidal debris, as the Jacobi radius in 
this snapshot is 36 pc. Hence the King template, which 
does not have a proper tidal debris term, gets significantly 
worse with increasing fit range. Furthermore, since there is 
an infiexion point of the slope between the bulk and the 
tidal debris at about the Jacobi radius, the templates which 
cannot account for the bulk and debris separately (EFF 
and Nuker) get worse for increasing fit range as well. In 
contrast to EFF, however, Nuker can partially compensate 
the missing bulk term by a more flexible core term. Besides 
Nuker, only the two templates which consist of three parts 
have a good reduced x^ value over all fit ranges (lower right 



panel of Fig. [5]). 

Taking a closer look at the tidal debris of this cluster 
at the given, dynamically well-evolved state shows that the 
slope outside the Jacobi radius is quite similar to the slope 
of the bulk, though it is not the same. If the data had been 
slightly more noisy we could have even concluded that the 
slope of the tidal debris close to th e cluster is the same as 
t he sl ope of the bulk, which is what lElson. Fall fc Freemanl 
(19871) might have seen in their data on young LMC 
clusters. This is also the reason why EFF and Nuker fit this 
cluster reasonably well even though they have no separate 
bulk and debris terms. 

Another striking fact is that the tidal debris does 
not show a constant slope for radii much larger than the 
Jacobi radius. From a constant mass-loss rate (which this 
cluster shows for several Gyr) and a constant velocity of 
the escapers within the tidal tails away from the cluster we 
would infer that the number of stars in the debris grows 
linearly with radial distance from the cluster. Thus, the 
surface density profile of the debris should fall off with R~^, 
as the surface area of the annuli grows with R^. Instead we 
see a slope which is steeper close to the cluster and gets 
more shallow at radii of several hundred pc. 

This is due to the fact that stars do not move linearly 
along the tidal tails as they escape from the cluster, but 
exhibit an epicyclic motion in which they get periodically 
accel erated and decelerated (|Kupper. Macleod fc Heggid 
[2003). Since escaping stars are first accelerated away from 
the cluster, the slope of the surface density close to the 
cluster should be lower than -1. This acceleration then 
turns into a deceleration at a distance R = 0.5 yc, where 
the length of the epicycles, yc, depends on the tidal field. 
At a distance R = yc, which in this case is given by ap- 
proximately SnRj ~ 340 pc, the slowest escapers (escapers 
with E — Ecrit — 0) are coming to a halt before they get 
re-accelerated. Escapers with a finite velocity above the 
escape velocity come to a halt at a somewhat larger yc, 
so that the mean escaper reaches its lowest drift velocity 
along the tidal tails at about 400-500 pc. (For a detailed 
description of the tidal tails of this particular cluster see 
' Kiipper et al.ll2010l ). In this range the slope of the surface 
density profile should be significantly larger than -1. 

Of greatest interest with regard to real star clusters is 
the tidal debris close to the clusters, as observations barely 
extend to radii of several Jacobi radii. Thus, we focus on 
the inner 100 pc in the following discussion. In Fig. [6] the 
evolution of the power-law slopes of the tidal debris can be 
seen as fitted by the given templates for a fit range of 1-100 
pc. We see that EFF and KKBH have a similar evolution, 
though KKBH varies less with time. This is due to the fact 
that EFF fits one slope for the bulk of the cluster and the 
tidal debris, whereas KKBH decouples the two parts and 
therefore gives a more accurate value for the slope of the 
tidal debris close to the cluster. As expected, the slope is 
significantly smaller than -1; indeed it has a value of about 
-5. Note that KKBH and EFF show a slightly shallower 
slope at the beginning of the simulation, where the cluster 
loses more mass as a consequence of primordial escapers 
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(see iKiipper et alj [2010|). By contrast, we see that the 
Nuker and the gen. Nuker templates do not give rehable 
information on the tidal debris slope. The power-law slopes 
of the Nuker and the gen. Nuker templates in Fig. [6] even 
had to be smoothed in order to reduce the stochastic noise 
to a reasonable level. Even so, the fluctuations in the values 
are too large to allow any conclusions to be drawn from 
them. This is due to the factors a and S in eq. [19] & [20] 
which are meant to adjust the smoothness of the transition 
from one slope to another but have quite a large influence 
on the fitting results. 

From Fig.[S]we can see that the surface density of bound 
stars falls off more steeply than the surface density of all 
stars in the sample (green profile vs. red profile), i.e. the 
edge of the bulk of bound stars lies further in than the edge 
of the bulk of all stars. Thus, the population of potential es- 
capers enhances or even dominates the surface density pro- 
file for radii larger than about half the Jacobi radius, just as 
we saw in the velocity dispersion profiles (Fig. [2]). ft should 
be kept in mind, however, that the underlying investigation 
is limited to star clusters with up to 64k stars. How this 
behaviour c hanges with increa sing A'^ still has to be consid- 
ered. From JBaumgardtl l|200ll '). however, we expect only a 
weak dependence on A'^. 

In Fig. [7] the time evolution of the Jacobi radius of this 
cluster is shown. As the cluster constantly loses mass the 
Jacobi radius (eq. llOf) slowly decreases. Also shown in the 
figure are the fitted edge radii of the King, the gen. Nuker 
and the KKBH templates for a fit range of 1-100 pc. We 
see that none of the templates reproduces the Jacobi radius 
accurately over the whole simulation time. 

The King Rt shows strong evolution at the beginning, 
because of an initial burst of mass-loss due to primordial 
escapers, which the template tries to fit. For the rest of the 
simulation the King edge radius agrees with the Jacobi ra- 
dius to within about 10%. 

The gen. Nuker Rt shows large fiuctuations, which are 
due to the above mentioned smoothness parameters a and S, 
which introduce a large ambiguity but also cause problems 
in the fitting process. Thus, for this template the Marquardt- 
Levenberg algorithm sometimes does not converge properly, 
and then the starting value for Rt, which was chosen to be 
90 pc for all templates, is retained. In case the fit converges, 
it mostly yields an edge radius which is about 50% too large. 

The KKBH template usually yields similar results to 
King, though 10% smaller, but is much less influenced by 
the initial mass loss of primordial escapers, and thus more 
accurately reproduces the edge radius of the bulk of the 
cluster. This case is a good example in which heavy mass 
loss signiflcantly influences fltting results, as may have 
happe ned in the investigation of JElson. Fall fc FreemanI 
lll987^■ or as was m entioned by ICote et all ( 20021 ). 
iCarraro. Zinn fc Moni Bidinl ((20o3) and lCarrarol (|2009l ) for 
the Milky- Way globular clusters Palomar 13, Whiting 1 and 
AM 4, respectively. The structure of the KKBH template 
helps to minimise such a problem. 

In Fig. [8] we show all these fltted edge radii in units 
of the Jacobi radius, averaged over all time steps of the 
simulation, for all 5 flt ranges. King shows a value which 
lies within 10% of the true Jacobi radius for large flt ranges. 



Only the smallest flt range of 1-25 pc yields an average edge 
radius which is about 30% off, so here the flt range would 
signiflcantly influence the fltting result. In contrast to that 
the KKBH template gives a more constant value over all 
flt ranges but which is reduced compared to King by about 
10%. The gen. Nuker template does not yield a reliable 
value for the edge radius, as the fluctuations between the 
time steps are large, as indicated by the error bars, and the 
mean strongly depends on the flt range. 

Throughout the analysis we found the bound stars to 
be well represented by a simple King template. From Fig. [5] 
we see that the edge of bound stars lies at about 70% of 
the Jacobi radius. The break radius of the KKBH template, 
at which it changes into the power-law tidal debris, is also 
shown in the flgure and lies at about 50% of Rj and can be 
interpreted as the radius at which the proflle of the cluster 
deviates from a smooth bulk proflle and gets dominated by 
potential escapers. 

In Tab. lBll a summary of the fltting results of the edge 
radii for the three smallest fltting ranges for all clusters 
in a constant tidal fleld is given. From this table we can 
conclude that the edge radii are closer to the Jacobi radius 
the larger the flt range is. For the smallest flt range the 
deviation can be even up to 40%. Averaging over all the 
fltting results in Tab. IBll conflrms our flndings: the King 
edge radius is on average a little (10-20%) smaller than the 
true Jacobi radius. Gen. Nuker yields no reliable results, 
whereas KKBH is more stable than King but yields an edge 
radius which is also about 20% smaller than the true Jacobi 
radius. Furthermore we see from the table that the edge of 
bound stars lies at about 70% of the Jacobi radius and that 
KKBH flnds a break in the proflle slope at about 50% of 
the Jacobi radius. Thus, at this radius the proflle becomes 
dominated by potential escapers. 

From the investigation of the clusters in a constant tidal 
fleld we can conclude that our ansatz for constructing a star 
cluster template consisting of core, bulk and debris was cho- 
sen well. The unbound tidal debris has a signiflcant influence 
on the fltting results, but this can be minimised by the ad- 
dition of an independent tidal debris term. In addition, our 
constructed KKBH template is a good tool with which to 
measure the power-law slope of the tidal debris and to look 
for a break in the surface density proflle. In this way we saw 
that the influence of potential escapers on the cluster proflle 
is signiflcant for radii larger than 0.5Rj and that beyond 
about 0.7Rj nearly all stars are unbound. 

4-2.2 Time- dependent tidal fields 

In Fig. [9] we see the clusters with an orbital eccentricity of 
0.25 (left panel) and 0.5 (right panel) at the time when they 
are at their apogalactic distance of 8.5 kpc, just as shown 
in Fig. [2] for the velocity dispersion proflles. Also shown in 
the flgure are the template flts for a flt range of 1-100 pc. 
We see no signiflcant difference from the surface density 
proflle of the cluster on a circular orbit (Fig. [5]), except for 
the mean density within the tidal debris, which gets larger 
with increasing eccentricity. This is due to the average mass 
loss of the clusters as well as the orbital compression of the 



Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters 15 



Data 

King 

EFF 

Nuker 

Gen. Nuker 

KKBH 

Data (bound stars) 

King (bound stars) 




1000 



R[pc] 



10'' 

10^ 

10^ 

cT 10° 

o 

Q. 

W 10"'' 



10"' 



10" 



10 



8.5 kpc 
£ = 0.5 




Data 
King 
EFF 
Nuker 
Gen. Nuker 
KKBH 
Data (bound stars) 
King (bound stars) 




1000 



R[pc] 



Figure 9. Surface density profiles of the clusters at an apogalactic distance of 8.5 kpc with an orbital eccentricity of 0.25 (Bl, left) and 
0.5 (B2, right). At the time of the snapshot both clusters are at apogalacticon and have a mass of about 15 000 Mq. The power-law slope 
of the tidal debris as measured by the KKBH template is »; = 5.2 it 0.3 on the left and rj = 4.5 ± 0.2 on the right. The inner cluster 
profiles look quite similar whereas the average density outside the cluster with an orbital eccentricity of 0.5 is higher, which is due to 
the higher mass-loss rate of this cluster and the larger orbital compression of its tidal tails. The theoretical Jacobi radii can be read off 
from the outermost green data points (bound stars). Corresponding velocity dispersion profiles are given in Fig. [21 
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Figure 11. Surface density profiles for the cluster at an apogalactic distance of 17 kpc with an orbital eccentricity of 0.75 (D3). On 
the left the cluster is shown shortly before reaching apogalacticon, in the middle it is at apogalacticon whereas on the right it is at 
perigalacticon. The power-law slope of the tidal debris as measured by the KKBH template is r; = 1.7 ± 0.2 on the left, r] = 2.9 ± 0.5 
in the middle and r] = 3.8 ± 0.1 on the right. Close to and at apogalacticon the cluster and its tidal tails get heavily compressed which 
leads to a great enhancement in the number of stars close to the cluster. At perigalacticon most of the profile is best fit by a power-law. 
Furthermore, most stars are unbound (green line), just as we have seen in the velocity dispersion profile (middle graph of Fig. [Sj. The 
drop in surface density in the left and middle panel for large radii is due to the young age of the cluster and its short tails, as escaped 
stars simply could not travel far at this time. Corresponding velocity dispersion profiles are given in Fig. |3] 



tails at apogalacticon, which both increase with eccentricity. 

The slope of the tidal debris close to the cluster 
(~ i?-") is rj = 5.2 ± 0.3 for e = 0.25 and r; = 4.5 ± 0.2 
for e = 0.5, i.e. about the same as in the constant tidal 
field (about t; = 5 for e — 0.0). But since the cluster has 
a lower orbital velocity at apogalacticon for increasing 
eccentricity, the epicyclic overdensities in the tidal tails are 
also located closer to the cluster (see iKiipper et al.l I2OI0I ) 
and the surface density is enhanced relative to the circular 



From Fig.[TO]we see that for a small orbital eccentricity 
of e = 0.25 the tidal debris slope evolves very much as 



in the constant tidal field case. For e = 0.5 we see that 
the slope on average shows the same evolution, but in 
addition shows jumps to smaller values. These jumps occur 
shortly before apogalacticon where the tails get maximally 
compressed due to the orbital acceleration, so that, for 
large eccentricities, they literally get pushed back inside the 
cluster. This enhances the tidal debris close to the cluster 
and thus temporarily yields a shallower slope. This effect 
is shown in Fig. [TT] for the cluster with an apogalactic 
distance of 17 kpc and an orbital eccentricity of 0.75. The 
left panel shows the cluster close to apogalacticon, the 
middle panel shows it at apogalacticon whereas the right 
panel shows the cluster at perigalacticon, similar to what 
we have seen in Fig. [31 The slopes rj of the three snapshots 
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Figure 10. Evolution of the slopes of the tidal debris of the 
templates as shown in Fig. \6\ for the clusters with an apogalactic 
distance of 8.5 kpc and an orbital eccentricity of 0.25 (Bl, upper 
panel) and 0.5 (B2, lower panel). The Nuker and the gen. Nuker 
slopes are smoothed with a cubic spline, as the variations between 
the individual time steps are very large. The EFF and KKBH 
(for e = 0.5 also the Nuker) slopes again show a similar evolution 
around an average value of about rj = 5. For increasing orbital 
eccentricity the power-law slope of the tidal debris starts to vary 
with the orbital phase. As the cluster gets compressed close to 
apogalacticon, the slope of the tidal debris abruptly decreases to 
values of about »; = 2 — 3. 



as measured by the KKBH template are rj = 1.7 ± 0.2 
close to apogalacticon, ri — 2.9 ± 0.5 at apogalacticon and 
7] = 3.8 ±0.1 at perigalacticon. 

Summarising the results, including also those from 
all the other models of our systematic study, which are 
not shown here in detail, we can conclude that the tidal 
debris close to the cluster (up to about three times the 
Jacobi radius) has a power-law slope ry in the range 4 to 
5, which can decrease to values in the range 1-2 shortly 
before apogalacticon due to orbital compression of the tidal 
tails. This may lead to the conclusion that we can identify 
those star clusters which are at an orbital phase close to 
apogalacticon by the slope of their surface density profiles 
outside the Jacobi radius. Star clusters like Palomar 13, 
AM 4 and Whiting 1 may well be in such an orbital phase, 
as their surf ac e den sit y profiles show slopes of about rj = 1.8 
dCote et al.] l2002l : ICarraro. Zinn fc Moni BidinI l2007l : 
ICarrardl2009l ). This assumption is further strengthened by 
the Milky- Way globular cluster Palomar 5, whose orbit is 
well constrained due to its long tidal tails (IDehnen et al.l 
[2OOJ). It is currently located close to apogalacti con and 
show s a power-law slope of about -1.5 IjOdenkirchen et al.l 
l2003l ). 



From the right panel of Fig. [TT] we can see that at 
perigalacticon, where the Jacobi radius is smallest (about 
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Figure 12. Evolution of the Jacobi radii, evaluated using eg. 1101 
(rod dotted line), and the edge radii of the fitted templates for a fit 
range of 1-100 pc for the clusters with an apogalactic distance of 
8.5 kpc and an orbital eccentricity of 0.25 (Bl, upper panel) and 
0.5 (B2, lower panel). The line colours are the same as in Fig. [7] 
Even though the Jacobi radius changes periodically the bulk of 
the cluster seems to be unaffected. The evolution of the King 
and KKBH edge radii (red and black solid line) are very much 
like in Fig. [7] but shifted to lower values for increasing orbital 
eccentricity. Only the break radius of the KKBH template (black 
dashed line) seems to be somewhat affected by the changing tidal 
field. 



10 pc) , most of the cluster is unbound (green data points vs. 
red data points) just as we have seen in the corresponding 
velocity dispersion profile (middle panel of Fig. |3]). Also 
note here how much more poorly the King template fits 
the cluster at perigalacticon than at apogalacticon. By 
contrast, the KKBH template, as well as the others with 
a power-law tidal debris term, can cover the whole cluster 
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and its debris perfectly well. 

In Fig. [12] we show the evolution of the Jacobi radius 
and of the fitted edge radii for the above clusters with 
eccentricities of 0.25 and 0.5, just as in Fig. [7] for e — 0.0. 
While the Jacobi radius oscillates with an amplitude which 
increases with increasing eccentricity, the overall evolution 
of the other radii is very similar to the constant tidal field 
case. The edge radius of the King template shows some 
evolution at the beginning but then gives the mean Jacobi 
radius of the cluster. The KKBH template again shows less 
evolution than King at the beginning, i.e. is less affected by 
heavy initial mass loss, but again gives a nearly constant 
value similar to that of King, but reduced by about 10%. 
The gen. Nuker edge radius does not show a well-behaved 
evolution at all. 

The results from all models on eccentric orbits in our 
set are presented in Tab. IB2IB4I sorted by eccentricity. The 
results of the clusters with e = 0.25 and e = 0.5 agree well 
with the findings for the clusters in a constant tidal field 
fTab. IBT)| . We note especially the consistent values for the 
edge of the bound stars at 70% of the Jacobi radius and 
the break radius of the KKBH template at 50% of the 
Jacobi radius. But also the fact that the King template 
yields on average about 90% of the mean Jacobi radius 
is interesting, since it is commonly thought that observed 
edge radii correspond to the perigalactic Jacobi radii of the 
corresponding clusters. The fact that the edge radius is not 
set at perigalacticon is due to the fact that stars do not 
leave the cluster immediately after becoming unbound but 
stay within the Jacobi radius for up to several orbits of the 
cluster about the galaxy. Moreover, when we consider the 
elongated shape of the equipotential surfaces of a cluster 
in a tidal field it is quite reasonable that a spherically 
symmetric, one-dimensional density profile will yield a 
smaller value than the two-dimensional (i.e. projected) 
Jacobi radius. 

Note that the mean results and fiuctuations get worse 
for increasing tidal field strength. All clusters with orbital 
eccentricities of 0.75 get disrupted within a few orbits about 
the galactic centre and thus are barely in a state close to 
virial equilibrium. Also the clusters at 4.25 kpc and 8.5 kpc 
with eccentricities of 0.5 show large spreads in the fits as 
they penetrate deeply into the central part of the galaxy. 
At the beginning of these A^-body computations, however, 
when the clusters are still very massive and much less 
vulnerable to tidal influences, the fitting results are more 
stable and agree with the findings for the other clusters. 



4-2.3 Concentrated clusters 

Fig. [13] shows the surface density profiles of the clusters at 
4.25 kpc (left) and 12.75 kpc (right), which both had an 
initial (3-dimensional) half-mass radius of 3 pc; for the first 
of these the corresponding velocity dispersion profile was 
shown in Fig. [3] At the time of the snapshots the clusters 
have 15 000 M0 and a ratio of projected half-mass radius 
to Jacobi radius of 0.08 (left) and 0.06 (right). From the 
figures and the fitted templates we see that the power-law 
slope of the tidal debris agrees with the previous findings 
(77 = 4.6 ± 0.6 for the left panel and ?? = 4.3 ± 0.3 for 
the right). Furthermore, we see that the EFF template 
completely fails to fit the data. This is due to the fact 
that the clusters do not have a distinct core (at least down 
to a radius of 1 pc) but appear to consist only of a bulk 
and tidal debris. Thus the EFF template, which consists 
of a flat core and a tidal debris term, cannot reproduce it. 
The King template does a better job for the bulk but fails 
at the edge of the bulk, whereas KKBH can give a good 
representation of all parts of both clusters. Nuker and gen. 
Nuker are also flexible enough to cover the whole radial 
range. These characteristics become more pronounced for 
increasing galactocentric radius as the half-mass radius is 
flxed at 3 pc but the Jacobi radius increases, and so the 
concentration increases. Note also that the quality of the 
King template flt decreases with a larger concentr ation. 
This supports the findings by Baumg ardt et al.l (|2010 ll that 
concentrated clusters cannot be properly fitted by a King 
template. 

But even though most templates yield good x^ values, 
from Tab. IB5I we can see that the fitted edge radii do not 
correspond to the Jacobi radii any more. With increasing 
concentration the edge radii shrink compared to the Jacobi 
radii. That is the reason why we avoid the term tidal radius 
and instead use edge radius, since this edge radius does not 
have to be a result of the galactic tide. 

Moreover, since the bound stars follow a similarly con- 
centrated distribution and since we tried fitting this distri- 
bution with a King template, these numbers also get unre- 
liable. Here a reconsideration of the templates with regard 
to the structure of concentrated clusters would be necessary, 
since the templates we used here all seem to fail in giving in- 
formation on the cluster structure, due to the fact that these 
clusters do not show a distinct core and bulk down to 1 pc. 
Thus, we conclude that Jacobi radii cannot be reliably deter- 
mined from observed clusters using these templates if their 
concentration is high, i.e. with Rhp/Rj — 0.1 or smaller. 



From the clusters in time-dependent tidal fields we can 
conclude that our KKBH template again proves to be a good 
tool to investigate the state of the clusters and their tidal 
debris. Thus, we were able to extend our findings of the 
previous section to clusters on eccentric orbits. Moreover, 
we suggest that star clusters in an orbital phase close to 
apogalacticon can be identified by the slope of their surface 
density profile outside the Jacobi radius. We furthermore 
find that the bulk of the cluster adjusts to the mean tidal 
field rather than, as usually believed, to the perigalactic tidal 
field. 



4.3 Review of templates 

Before summarising the results of this investigation let us 
briefly discuss the performance of the templates which were 
defined in Sec. 13.41 and which were used here. 



• King has proved to be a reasonable tool for studying 
the Jacobi radius of extended clusters but fails when the 
cluster has a pronounced tidal debris component or is too 
concentrated. It cannot give information on the core slope 
of either mass segregated clusters or clusters with an IMBH 
in the centre. Eccentric cluster orbits are also a problem for 
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Figure 13. Surface density profiles of the concentrated clusters at a galactocentric radius of 4.25 kpc (AOc, left) and 12.75 kpc (COc, 
right). At the time of the snapshots both clusters have a mass of about 15 000Mq. The power-law slope of the tidal debris as measured 
by the KKBH template is rj = 4.6 ± 0.6 (left) and t] = 4.3 ± 0.3 (right). The theoretical Jacobi radii can be read off from the outermost 
green data points (bound stars). One corresponding velocity dispersion profile is given in Fig. |4] 



the template. Moreover, the fit results strongly depend on 
the radial range of the cluster which is covered; it is thus 
quite unreliable for obtaining cluster parameters. 

• EFF is a good tool for the study of clusters with a 
pronounced tidal debris component, provided that they are 
not too concentrated or mass segregated. It cannot, however, 
give information on the edge radius. 

• Nuker can fit most of the clusters in our sample well, 
but it is not possible to deduce cluster parameters reliably 
from the fitting results. Furthermore it does not have an 
edge radius. 

• Gen. Nuker is a very powerful template with a low 
X^ for all fit ranges but, just as Nuker, does not allow the 
deduction of reliable structural parameters from the fitting 
results, as the resulting parameters of the template vary 
strongly and depend strongly on the radial coverage of the 
cluster. 

• KKBH has proved to be a reasonable extension of 
King, improving the fitting results for all clusters in the 
sample. It also yields the most reliable structural informa- 
tion over all fit ranges in terms of edge radius and the slope 
of the tidal debris. Like all other templates, however, it fails 
in the case of concentrated clusters. 



5 CONCLUSIONS 

This systematic investigation of velocity dispersion profiles 
and surface density profiles of star clusters with masses of 
about 10* Mq has shown how complex and versatile these 
profiles can be, depending on the circumstances of the 
clusters. 

Furthermore we have shown the importance of po- 
tential escapers for these two quantities. Both profiles are 
influenced by this population of energetically unbound stars 
for radii larger than about 50% of the Jacobi radius and 



are completely dominated by it for radi i larger than about 
70% of the Jacobi radius (Fig. [2] & [Sjl. iBaumgardtl (|200ll ) 
found that for clusters in a constant tidal field the fraction 
of potential escapers varies with the number of stars within 
a cluster approximately as A'^"^'*. Thus, even if the cluster 
mass increases by a factor of 10, the fraction (and infiuence) 
of potential escapers will only be reduced by less than 50% 
compared to this investigation. 

As a consequence of the potential escaper population, 
the velocity dispersion profiles in our sample do not show 
a clear separation between the bulk of the cluster and the 
tidal debris at the Jacobi radius, but rather show a smooth 
transition from cluster to debris. For mass estimates which 
are based on velocity dispersion measurements we therefore 
recommend to use stars within 50% of the Jacobi radius to 
minimise the effect of potential escapers. 



Moreover, we 
locity dispersion 



argue that investigations onve- 

profiles li ke iDrukier et all ll 19981) 

Scarpa. Marconi fc Gilmozzil (|2003l ) or I Scarpa et al.l 



axp, 

m 



20071 ) have to be interpreted with caution. Also, detecting 



a possible dark matter halo around globular clusters is less 
feasible due to the population of potential escapers. Based 
on our investigation we suggest that no deviation from 
Newtonian gravity in a star cluster has been detected so 
far, neither has there been evidence for a dark matter halo. 

From our large set of surface density profiles we saw 
that in most cases the structure of extended star clusters 
can be split into core, bulk and tidal debris. 

(i) The core is the innermost part of the cluster which 
extends out to the core radius. In most cases it is fiat, 
and only for clusters with strong mass segregation or with 
an IMBH is it cuspy. Very concentrated clusters with 
Rhp/Rj < 0.1 have such a small core that fitting such 
clusters leads to problems with templates as well as with 
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physical models (see also iBaumgardt et al.ll201(]| ). 



(ii) The bulk contains most stars of the cluster and 
extends hom the core to the tidal debris. We found that 
the bulk can be well represented by a King-like template. 
We named the radius at which such a template approaches 
zero density the edge of the bulk and not the tidal radius of 
the cluster, since we found it not to be consistent with the 
Jacobi radius. In the case of concentrated clusters, the edge 
radius is not even determined by tidal forces at all fTab. lBS)) . 

Moreover, for clusters on eccentric orbits we found that 
the edge radius of a cluster adjusts to the mean Jacobi 
radius (Fig. I12p and not the perigalactic Jacobi radius, 
which has been assumed in earlier investigations (e.g. 
Ilnnanen. Harris, fc Wobbink 1983, Fall & Zhang 2001), but 
has never been checked by self-consistent calculations on a 
star-by-star basis. Since nearly all globular clusters of the 
Milky Way are more compact and massive than our test 
cluster, and thus should be less influenced by tidal varia- 
tions, this finding should also hold for most of these clusters. 

The edge radius of extended and moderately concen- 
trated clusters, as fitted by the King template, lies within 
about 10% accuracy at 90% of the mean Jacobi radius 
(Tab. [BTIlB4p . But we found the King template to be a bad 
representation of the true profile in the case of non-circular 
cluster orbits because it is significantly influenced by 
the tidal debris and by the potential escaper population. 
Furthermore the accuracy of the results strongly depends 
on the radial coverage of the cluster. We therefore created 
an enhanced version of the King template which has more 
flexibility in the core and has an additional tidal debris 
term, given by a power-law. With this KKBH template 
we achieved a higher stability over all fit ranges. Further- 
more, we were able to measure the edge of the bulk more 
accurately and found it to lie at 80% of the Jacobi radius. 
This discrepancy is easily understandable if we take into 
account the fact that the edge radius of a template fits 
the azimuthally averaged mean of a cluster's tidal surface, 
whereas the Jacobi radius is by definition the semi-major 
axis of the equipotential surface. 



suggest that this enhanced slope can be used to identify 
star clusters which are close to or at the apogalacticon of 
their orbit. The most prominent example of such a cluster 
is the Milky Way globular cluster Palomar 5 which is 
well known to be c lose to the apogalacticon of its orbit 
(|Dehnen et al.ll2004r ). and which shows a power-law slope of 
-1.5 outside its assumed Jacobi radius (|Odenkirchen et al.l 
|2003). Moreover, the MW globular clusters Palomar 13, AM 
4 and Whiting 1 may be further candidates for being close 
to apogalacticon, since their surface density profiles show 
power-law slopes of about -1.8 outside the Jacobi radius, 
which has been attributed to heavy ongoing m ass loss by 
the authors of the corresponding i nvestigations (C6teet_aL| 
l2002l : ICarraro. Zinn fc Moni Bidinl.2007. : .Carraro..2009l )^ 

Furthermore we showed that due to the population of po- 
tential escapers the power-law slope of the tidal debris be- 
gins even before the Jacobi radius is reached. Thus, the bulk 
and the tidal debris partially overlap. With our KKBH tem- 
plate we found the debris to dominate for radii larger than 
about 50% of the Jacobi radius (Tab. lBUBl)) . 

With the help of our comprehensive set of clusters it 
is possible to tailor a star cluster template to the facts we 
have found. Our KKBH template seems to be a good first 
step but, like all other available templates, shows discrep- 
ancies for concentrated clusters. Of course, templates are 
less attractive than physically motivated models. None of 
the existing models, like King (1966) and Wilson (19751), 
however, account for the above mentioned facts. We suggest 
that a general, physically motivated model for star clusters 
should include a term for mass segregation and should 
include a potential escaper population when being fitted to 
observations. Finally, the possible influence of (compressed) 
tidal tails on the surface density proflles of star clusters 
should be kept in mind when fltting surface density profiles. 

Moreover, the fact that all clusters in the sample showed 
a constant edge of the bulk with time (Fig. I12|l and that 
stars which get unbound escape from a cluster with a delay, 
supports our theoretical treatment of epicyclic overdensities 
in tidal tails f or clu sters in time-dependent tidal fields in 
iKiipper et"aLll|2O10D . 



Furthermore we found that the edge radii from concen- 
trated clusters with Rhp/Rj < 0.1 cannot be extracted 
properly with currently available templates. The higher 
the concentration of the cluster the smaller is the ratio of 
the fitted edge radius to the theoretically evaluated Jacobi 
radius (Tab. lBS)) . 
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(iii) The tidal debris falls off like a power-law with a 
slope of about -4 to -5, rather than the expected dependence 
of R~^. This is due to the epicyclic motion of the stars in 
the tidal tails, which is most pronounced in the vicinity of 
the cluster. 

For clusters on eccentric orbits, at a time shortly before 
reaching apogalacticon, the slope was found to deviate 
for a short time to a value of about -2 as a consequence 
of orbital compression of the tails (Fig. [10] & lllf) . We 
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APPENDIX A: FITTING KKBH USING 
GNUPLOT 

When using GNUPLOT, fitting the KKBH template is quite 
simple since it is possible to use ternary operators. First 
the two functions fi{R) and f2{R) feQ.I22l& l23|) need to be 
defined with 

fl(x) = (x<Rt ? k * (x/Rc/(1.0+x/Rc))**(-gaimna) 
* (1.0/(1.0+(x/Rc)**2.0)**(0.5) - 
1.0/(1.0+(Rt/Rc)**2.0)**(0.5))**2.0 : 0) 

and 

f2(x) = fl(Rt*mu) * (1.0 + 

(x/ (Rt*mu) ) **64 . 0) ** (-eta/64 . 0) 

after which the KKBH function is defined 

f(x) = (x < mu*Rt ? fl(x) + b : f2(x) + b) 

is defined. Then the constants have to be set to some 
reasonable value, e.g. 

k = 1.0 
Re = 5.0 
gamma = 0.01 
Rt = 50.0 
mu = 0.5 
eta =3.0 
b = 0.0001 

If surface density data exists in a file named, e.g., 
data.txt, with columns R, S(_R) and AE(_R), the final 
fitting is carried out by use of the command 

fit f(x) /data.txt/ u 1:2:3 via k, Re, gamma, Rt , 
eta, mu, b. 

APPENDIX B: DATA TABLES 



This paper has been typeset from a T^^X/ WF^ file prepared 
by the author. 
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Table Bl. Fitting results for the edge radii in units of the theoretical Jacobi radius for the clusters in a constant tidal field 
averaged over all available snapshots of the specific model. The results are given for the three smallest fit ranges. In the last line 
a weighted mean and the standard deviation of the results is given. Rt gives the edge radius of the specific template, whereas Rj 
stands for the theoretical Jacobi radius, fi is the fraction of the edge radius, at which the KKBH template turns into a power-law. 



Name 


fit range 


Kmg 


Gen. Nuker 


KKBH 


King (bound) 




[pc] 


Rt 
Rj 


Rt 
Rj 


Rj 


Rt 
Rj 


Rt 
Rj 


AO 


1-25 


0.86 ± 0.06 


3.70 ± 0.86 


0.49 ± 0.09 


0.82 ± 0.07 


0.77 ± 0.47 




1-50 


0.97 ± 0.06 


2.12 ± 0.98 


0.53 ± 0.08 


0.85 ± 0.07 


0.74 ± 0.13 




1-100 


0.98 ± 0.06 


2.22 ± 0.97 


0.57 ± 0.11 


0.86 ± 0.08 


0.74 ± 0.13 


BO 


1-25 


0.68 ± 0.05 


2.58 ± 0.41 


0.41 ± 0.12 


0.75 ± 0.06 


0.70 ± 0.04 




1-50 


0.89 ± 0.05 


1.48 ± 0.70 


0.49 ± 0.06 


0.81 ± 0.06 


0.71 ± 0.04 




1-100 


0.92 ± 0.07 


1.36 ± 0.49 


0.51 ± 0.07 


0.81 ± 0.06 


0.71 ± 0.04 


CO 


1-25 


0.63 ± 0.04 


1.94 ± 0.35 


0.60 ± 0.05 


0.76 ± 0.03 


0.67 ± 0.02 




1-50 


0.84 ± 0.01 


1.56 ± 0.89 


0.46 ± 0.05 


0.79 ± 0.03 


0.71 ± 0.01 




1-100 


0.92 ± 0.05 


1.36 ± 0.74 


0.50 ± 0.05 


0.81 ± 0.03 


0.71 ± 0.01 


DO 


1-25 


0.56 ± 0.05 


1.58 ± 0.26 


0.57 ± 0.02 


0.71 ± 0.02 


0.63 ± 0.01 




1-50 


0.78 ± 0.01 


1.32 ± 0.96 


0.39 ± 0.03 


0.76 ± 0.02 


0.71 ± 0.01 




1-100 


0.90 ± 0.05 


1.20 ± 0.74 


0.49 ± 0.05 


0.80 ± 0.03 


0.71 ± 0.03 


mean 




0.81 ± 0.14 


1.84 ± 0.72 


0.50 ± 0.06 


0.78 ± 0.15 


0.69 ± 0.04 



Table B2. Fitting results for the edge radii in units of the theoretical Jacobi radius for the clusters in time-dependent tidal fields 
and orbital eccentricities of 0.25 averaged over all available snapshots of the specific model. The results are given for the three 
smallest fit ranges. In the last line a weighted mean and the standard deviation of the results is given. Rt gives the edge radius 
of the specific template, whereas Rj stands for the theoretical Jacobi radius. ^ is the fraction of the edge radius, at which the 
KKBH template turns into a power-law. 



Name 


fit range 


King 


Gen. Nuker 


KKBH 


King (bound) 




[pc] 


Rt 
Rj 


-Rt 
Rj 


t^Rt 
Rj 


Rt 
Rj 


Rt 
Rj 


Al 


1-25 


0.90 ± 0.19 


3.83 ± 1.76 


0.56 ± 0.13 


0.84 ± 0.17 


0.71 ± 0.42 




1-50 


0.98 ± 0.22 


3.37 ± 1.51 


0.58 ± 0.13 


0.86 ± 0.17 


0.69 ± 0.04 




1-100 


0.98 ± 0.22 


2.08 ± 1.59 


0.65 ± 0.14 


0.88 ± 0.18 


0.69 ± 0.04 


Bl 


1-25 


0.76 ± 0.12 


2.97 ± 0.70 


0.44 ± 0.11 


0.80 ± 0.13 


0.71 ± 0.04 




1-50 


0.95 ± 0.16 


1.60 ± 0.69 


0.54 ± 0.10 


0.86 ± 0.14 


0.72 ± 0.08 




1-100 


0.97 ± 0.17 


1.54 ± 0.57 


0.55 ± 0.12 


0.85 ± 0.15 


0.72 ± 0.08 


CI 


1-25 


0.66 ± 0.06 


2.18 ± 0.41 


0.63 ± 0.09 


0.81 ± 0.08 


0.68 ± 0.03 




1-50 


0.90 ± 0.09 


1.59 ± 0.94 


0.51 ± 0.07 


0.85 ± 0.08 


0.72 ± 0.02 




1-100 


0.96 ± 0.11 


1.38 ± 0.53 


0.54 ± 0.07 


0.86 ± 0.08 


0.72 ± 0.02 


Dl 


1-25 


0.69 ± 0.13 


1.86 ± 0.49 


0.63 ± 0.11 


0.79 ± 0.13 


0.66 ± 0.05 




1-50 


0.87 ± 0.15 


1.49 ± 1.00 


0.47 ± 0.09 


0.82 ± 0.14 


0.71 ± 0.03 




1-100 


0.97 ± 0.18 


1.30 ± 0.53 


0.53 ± 0.09 


0.85 ± 0.14 


0.71 ± 0.03 


mean 




0.85 ± 0.12 


1.92 ± 0.84 


0.55 ± 0.06 


0.84 ± 0.03 


0.70 ± 0.02 
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Table B3. Fitting results for the edge radii in units of the theoretical Jacobi radius for the clusters in time-dependent tidal fields 
and orbital eccentricities of 0.5 averaged over all available snapshots of the specific model. The results are given for the three 
smallest fit ranges. In the last line a weighted mean and the standard deviation of the results is given. Rt gives the edge radius 
of the specific template, whereas Rj stands for the theoretical Jacobi radius, fi is the fraction of the edge radius, at which the 
KKBH template turns into a power-law. 



Name 


fit range 


King 


Gen. Nukcr 


KKBH 


King (bound) 




[pc] 


Rt 
Rj 


-Rt 
Rj 


Rj 


-Rt 
Rj 


-Rf 
Rj 


A2 


1-25 


1.01 ± 0.53 


3.39 ± 2.75 


0.54 ± 0.25 


0.83 ± 0.40 


0.61 ± 0.10 




1-50 


1.18 ± 0.68 


2.75 ± 2.45 


0.58 ± 0.27 


0.95 ± 0.68 


0.62 ± 0.10 




1-100 


1.23 ± 0.73 


3.46 ± 2.49 


0.68 ± 0.56 


1.08 ± 1.04 


0.62 ± 0.10 


B2 


1-25 


0.80 ± 0.31 


2.72 ± 1.72 


0.50 ± 0.19 


0.79 ± 0.30 


0.66 ± 0.07 




1-50 


0.94 ± 0.37 


2.30 ± 1.63 


0.51 ± 0.18 


0.81 ± 0.31 


0.66 ± 0.07 




1-100 


0.97 ± 0.38 


1.58 ± 1.29 


0.55 ± 0.19 


0.82 ± 0.31 


0.66 ± 0.07 


C2 


1-25 


0.70 ± 0.24 


2.59 ± 1.18 


0.45 ± 0.17 


0.76 ± 0.25 


0.65 ± 0.11 




1-50 


0.90 ± 0.32 


1.92 ± 0.98 


0.51 ± 0.17 


0.82 ± 0.28 


0.67 ± 0.07 




1-100 


0.95 ± 0.34 


1.59 ± 1.02 


0.54 ± 0.18 


0.82 ± 0.28 


0.67 ± 0.07 


D2 


1-25 


0.64 ± 0.22 


2.25 ± 0.88 


0.52 ± 0.19 


0.74 ± 0.25 


0.62 ± 0.10 




1-50 


0.87 ± 0.31 


1.57 ± 0.91 


0.49 ± 0.16 


0.79 ± 0.27 


0.66 ± 0.07 




1-100 


0.94 ± 0.35 


1.33 ± 0.74 


0.55 ± 0.18 


0.81 ± 0.27 


0.66 ± 0.07 


mean 




0.88 ± 0.17 


2.06 ± 0.72 


0.52 ± 0.06 


0.81 ± 0.09 


0.65 ± 0.02 



Table B4. Fitting results for the edge radii in units of the theoretical Jacobi radius for the clusters in time-dependent tidal fields 
and orbital eccentricities of 0.75 averaged over all available snapshots of the specific model. The results are given for the three 
smallest fit ranges. In the last line a weighted mean and the standard deviation of the results is given. Rt gives the edge radius 
of the specific template, whereas Rj stands for the theoretical Jacobi radius, fi is the fraction of the edge radius, at which the 
KKBH template turns into a power-law. 



Name 


fit range 


King 


Gen. Nuker 


KKBH 


King (bound) 




[pc] 


-Rt 
Rj 


Rt 
Rj 


I^Rt 
Rj 


Rt 
Rj 


-Rf 
Rj 


A3 


1-25 


2.06 ± 2.01 


4.99 ± 2.44 


1.16 ± 2.23 


1.08 ± 0.72 


1.47 ± 2.06 




1-50 


2.26 ± 2.27 


3.67 ± 2.31 


1.20 ± 1.47 


2.07 ± 2.44 


1.43 ± 1.90 




1-100 


2.17 ± 2.09 


4.12 ± 2.05 


1.91 ± 2.23 


1.31 ± 1.28 


1.43 ± 1.90 


B3 


1-25 


0.64 ± 0.22 


2.25 ± 0.88 


0.52 ± 0.19 


0.74 ± 0.25 


0.62 ± 0.10 




1-50 


0.87 ± 0.31 


1.57 ± 0.91 


0.49 ± 0.16 


0.79 ± 0.27 


0.66 ± 0.07 




1-100 


0.94 ± 0.35 


1.33 ± 0.74 


0.55 ± 0.18 


0.81 ± 0.27 


0.66 ± 0.07 


C3 


1-25 


1.38 ± 1.39 


3.12 ± 2.11 


0.26 ± 1.68 


1.60 ± 1.78 


0.95 ± 1.21 




1-50 


1.96 ± 1.81 


2.57 ± 1.94 


0.88 ± 1.06 


1.57 ± 1.53 


0.85 ± 0.73 




1-100 


2.24 ± 1.83 


2.60 ± 1.91 


1.04 ± 1.40 


1.77 ± 1.84 


0.87 ± 0.76 


D3 


1-25 


1.36 ± 1.67 


2.67 ± 1.85 


1.17 ± 1.73 


1.43 ± 1.53 


1.07 ± 1.15 




1-50 


1.75 ± 1.66 


2.09 ± 1.73 


1.25 ± 1.94 


1.74 ± 1.89 


1.25 ± 1.47 




1-100 


2.06 ± 1.74 


2.10 ± 1.64 


1.36 ± 1.80 


1.74 ± 1.65 


1.29 ± 1.52 


mean 




1.14 ± 0.58 


2.43 ± 1.06 


0.67 ± 0.46 


1.03 ± 0.44 


0.72 ± 0.32 
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Table B5. Fitting results for the edge radii in units of the theoretical Jacobi radius for the concentrated clusters with a half-mass 
radius of 3 pc averaged over all available snapshots of the specific model. The results are given for the three smallest fit ranges. In 
the last line a weighted mean and the standard deviation of the results is given. Rt gives the edge radius of the specific template, 
whereas Rj stands for the theoretical Jacobi radius, /i is the fraction of the edge radius, at which the KKBH template turns into 
a power-law. 



Name 


fit range 


King 


Gen. Nukcr 


KKBH 


King (bound) 




[pc] 


Rt 
Rj 


-Rt 
Rj 


M-Rt 

Rj 


-Rt 
Rj 


Rt 
Rj 


AOc 


1-25 


0.75 ± 0.10 


3.72 ± 0.24 


0.35 ± 0.10 


0.65 ± 0.11 


0.79 ± 0.92 




1-50 


0.79 ± 0.10 


1.85 ± 0.93 


0.43 ± 0.10 


0.72 ± 0.11 


0.69 ± 0.06 




1-100 


0.80 ± 0.10 


2.46 ± 0.90 


0.44 ± 0.12 


0.71 ± 0.12 


0.69 ± 0.06 


BOc 


1-25 


0.56 ± 0.10 


2.51 ± 0.25 


0.36 ± 0.23 


0.55 ± 0.19 


0.59 ± 0.10 




1-50 


0.64 ± 0.12 


2.17 ± 0.90 


0.33 ± 0.13 


0.57 ± 0.14 


0.64 ± 0.43 




1-100 


0.64 ± 0.13 


1.50 ± 0.51 


0.33 ± 0.09 


0.58 ± 0.14 


0.64 ± 0.43 


COc 


1-25 


0.47 ± 0.12 


1.87 ± 0.60 


0.34 ± 0.19 


0.49 ± 0.19 


0.52 ± 0.12 




1-50 


0.55 ± 0.14 


1.48 ± 0.57 


0.20 ± 0.06 


0.43 ± 0.12 


0.53 ± 0.13 




1-100 


0.56 ± 0.15 


1.39 ± 0.34 


0.26 ± 0.10 


0.49 ± 0.15 


0.53 ± 0.13 


DOc 


1-25 


0.36 ± 0.06 


1.65 ± 1.04 


0.26 ± 0.15 


0.38 ± 0.14 


0.41 ± 0.07 




1-50 


0.42 ± 0.09 


1.36 ± 0.81 


0.14 ± 0.02 


0.32 ± 0.06 


0.42 ± 0.09 




1-100 


0.43 ± 0.10 


1.29 ± 0.32 


0.17 ± 0.05 


0.36 ± 0.08 


0.42 ± 0.09 


mean 




0.57 ± 0.15 


2.10 ± 0.70 


0.24 ± 0.10 


0.50 ± 0.13 


0.55 ± 0.12 



